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Abstract. It is assumed that a black hole can be disturbed in such a way that a 
ringdown gravitational wave would be generated. This ringdown waveform is well 
understood and is modelled as an exponentially damped sinusoid. In this work we use 
this kind of waveform to study the performance of the SCHENBERG gravitational 
wave detector. This first realistic simulation will help us to develop strategies for the 
signal analysis of this Brazilian detector. We calculated the signal-to-noise ratio as 
a function of frequency for the simulated signals and obtained results that show that 
SCHENBERG is expected to be sensitive enough to detect this kind of signal up to a 
distance of ~ 20kpc. 
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1. Black hole ringdown 



Ringdown waveforms originate from a small perturbation of a spinning black hole (BH) 
and can be modelled as an exponential-damped sinusoid. The central frequency / of 
the fundamental quadrupolar mode and the quality factor Q depend on the mass M 
and the spin (S = aGM 2 /c) of the BH. They can be approximated by an analytic fit 
[U |2] and are given by 



/ « 32kHz 



i -0.63(1 - a; 



M 



M 



and Q ~ 2 (1 
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Figure 1. Volume formed by sets of parameters that correspond to frequencies varying 
from / = 3172.5Hz to 3240Hz. 



The value of the dimensionless spin parameter a is zero in the Schwarzschild limit 
(non- rotating BH) and one in the extreme Kerr limit (maximum rotational speed), so 
Q > 2. 

The relative strain caused by any metric perturbation on the detector depends 
on the position of the BH in the sky and the relative orientation of its spin axis to 
the local zenith of the detector. The averaged strain produced by the tensorial GW 
components can be obtained by rms averaging over all possible angles and the result is 
(^Gw)angie = Aq(t) , where q{t) is the damping function that depends on the rotation 
speed and on the central frequency of the BH: 

q(t) = {27r)h- nft/Q cos{27c ft) (2) 

for t > (we set the time origin on t = 0). The amplitude A also depends on the 
fraction e of the total mass-energy radiated and on the BH distance from the Earth (r), 
and it is given by 

Mpc\ f_e 

The quality factor roughly gives the number of coherent cycles for the waveform 
It affects the damping (1/Q), which decreases with increasing rotational speed. 
To generate the signal that was introduced in the model of the detector we chose one 
(among many) of the sets of parameters that made the central frequency of the BH 
quadrupolar normal mode coincident with one of SCHENBERG resonant quadrupolar 
frequencies (3172. 5. 6Hz, 3206.3Hz and 3240Hz). We obtain a certain volume in the 
space of parameters shown in figure Q by fixing two boundary frequencies for the BH 
quadrupolar mode and then making those parameters to vary. Consequently any set of 
parameters inside that volume is within SCHENBERG bandwidth. 



A ~ 2.415xlO" 21 Q- ; 



1-0.63(1 



a) w 



.(3) 
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Figure 2. Waveform of a black hole ringdown (a) in the time domain and (b) in the 
frequency domain. 

2. The simulated signal 

We chose e = 0.01, M = 4M and a = 0.16 to generate the signal that we used to 
calculate the signal-to-noise ratio. We believe that this efficiency is not overestimated 
since similar values have been found in the literature [31 H] . 

Figure EJa shows the waveform h(t) in the time domain. The Fourier transform of 
the impulse h(t) is defined as 

/oo 
h(t)e- iuJt dt. (4) 
-oo 

It represents the signal in the frequency domain and is shown in Figure |21b. This figure 
shows that the signal spectrum on SCHENBERG bandwidth (3.17 — 3.24kHz) is almost 
plane. This is a direct consequence of the low value for the BH's quality factor (Q « 2.16, 
for that simulated signal). For higher values of the quality factor the damping function 
tends to a sinusoid and a peak may appear in the spectrum. This would facilitate the 
detection of the signal. 

3. The SCHENBERG model 

We modelled the SCHENBERG detector by assuming a linear elastic theory. Adopting 
this approach we determined the mechanical response of the system and obtained an 
expression for the case when six 2-mode mechanical resonators are coupled to the 
antenna's surface according to the arrangement suggested by Johnson and Merkowitz, 
the truncated icosahedron configuration E|. We found an analytic expression to 
calculate the spectral density for the mode channels, given by [Zj 

W = IU")| a sf + £ \a*(u)\%"* + J2 \a%(<»)\%"*- (P) 

5=1 j=l 

The scalar £ m (u;) represents the transfer function of the sphere's internal forces 
at angular frequency oo on mode m. The matrix elements VL^Auj) correspond to the 
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Figure 3. Sensitivity curve of SCHENBERG at 4.2K for the quadrupolar normal 
modes m. 



response functions of channel m to the noise forces F- ' that are acting on the mode 
i of resonator j. S x denotes the spectral density of the quantity X. In order to find 
equation EJ we assumed that all noise sources and possible signals were statistically 
independent so the cross terms could be neglected. 

The GW effective force El E] that acts on the sphere (in the frequency domain) 
has spectral density 

2 

(6) 



St GW (u) 



l -u 2 m sX B) ,v;;,u<). 



The sensitivity curve for mode m is given through 

h m {u) = 



1 



C9 



12' 



(7) 



{\uj 2 m s xR) 2 l&«MI 

Figure El shows the sensitivity curves for the five quadrupolar normal modes of the 
sphere. We obtained this curve assuming that no signal was present and that only 
the usual noises in this kind of detector (Brownian noise, serial noise, back-action 
noise and electronic phase and amplitude noises) \ were disturbing the antenna at 
4.2K. The absence of a signal made filtering unnecessary in this calculation. Also, 
we assumed a perfectly symmetric sphere, which implies identical transfer functions 
£ m (u;). Asymmetries on the sphere could make the transfer function of a mode different 
from the transfer function of the others. 

X For details on the parameters used in the Schenberg model see Aguiar, O. D. et al. and Frajuca, C. 
et al. in these Proceedings. 
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4. Signal-to-noise ratio 

Using the results in sections El (which inform about the signal) and El (which inform 
about the noise) it is easy to calculate the integrated signal-to-noise ratio (SNR). From 
the literature [10] we have 



The integral in equation |H1 which is dominated by the interval of SCHENBERG 
bandwidth, gives SNR > 1 for the simulated signal when the BH's distance is r ~ 20kpc. 

5. Conclusion 

If a low-spinning BH (a ~ 0.16) with mass ~ 4M and emission of about 1% of its mass- 
energy in the form of GW exists up to a distance of ~ 20kpc, then the SCHENBERG 
detector is expected to observe it when operating at 4.2K. 

In the calculation we didn't take into account other possible dependencies on the 
BH parameters except the known frequency-mass dependence. Strong dependencies 
between e and a may make signals from high-spinning BHs much stronger. 

There are theories that claim the existence of BHs but experimental evidences 
are indirect and inconclusive. Certainly it is not yet possible to predict the correct 
number of galactic BHs. However, the radius of ~ 20kpc includes the galactic center 
and the majority of the spiral arms, which are the most probable birthplaces for stellar 
BHs. So the detection of the kind of signal investigated in this work is an opportunity 
for gravitational wave astronomers to make essential discoveries about the black hole 
population in our galaxy. 
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